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Inclusive dijet production at large pseudorapidity intervals (An) between the two jets has been 
suggested as a regime for observing BFKL dynamics. We have measured the dijet cross section for 
large Ar; in pp collisions at = 1800 and 630 GeV using the D0 detector. The partonic cross 
section increases strongly with the size of An. The observed growth is even stronger than expected 
on the basis of BFKL resummation in the leading logarithmic approximation. The growth of the 
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partonic cross section can be accommodated with an effective BFKL intercept of Qbfkl (20 GeV) = 
1.65 ±0.07. 
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Jet production in the high-energy limit of Quantum 
Chromodynamics (QCD), as defined by center-of-mass 
(cm.) energies (\/~s) much larger than the momentum 
transfers (Q), presents a very interesting and yet little 
explored area. In this kinematic region, the significantly 
different energy scales of the process lead to calculat- 
ed jet cross sections characterized by the appearance of 
large logarithms ln(s/Q 2 ), which must be summed to all 
orders in a s . This summation is accomplished through 
the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [jlj equation, 
which involves a space-like chain of an infinite number of 
gluon emissions. The gluons have similar transverse mo- 
menta, but they are strongly ordered in their pseudora- 
pidities or, equivalently, in their longitudinal momentum 
fractions, x*. Thus, the BFKL equation effectively de- 
scribes the evolution in x (growth with 1/x) of the gluon 
momentum distribution in the proton. 

Attempts to isolate and probe the BFKL evolution in 
the low- a; region in ep collisions at HERA, by measuring 
the forward jet and particle cross sections have led 
to ambiguous results. (At HERA, forward denotes the 
region away from the current jet and towards the proton 
remnant.) In pp collisions, the azimuthal decorrelation 
as a function of the pseudorapidity interval, A77, in dijet 
systems has been studied ||. (Here, 77 = — ln[tan(0/2)], 
where 6 is the polar angle of the jet relative to the pro- 
ton beam.) It has been argued however, that the 
azimuthal angle distribution is not an inclusive enough 
quantity for the cancellation between the real and virtual 
soft gluons that accompany dijet production; thus, such 
a distribution can not be used as a probe of the BFKL 
equation. 

Inclusive dijet production at large pseudorapidity in- 
tervals in high energy pp collisions, on the other hand, 
provides an excellent testing ground for BFKL dynam- 
ics. We present a measurement of the dijet cross section 
at large A77 using the D0 detector at the Fermilab Teva- 
tron collider. We reconstruct the event kinematics using 
the most forward/backward jets, and measure the cross 
section as a function of x%, X2 and Q 2 . The longitudi- 
nal momentum fractions of the proton and antiproton, X\ 
and X2 , carried by the two interacting partons are defined 
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where Et x (Et 2 ) and 771 (772) are the transverse energy and 
pseudorapidity of the most forward(backward) jet, A77 = 
?7i — f/2 > 0, and fj — (rji+rfe) /2. The momentum transfer 
during the hard scattering is defined as: 



Q — v E Tl E T2 
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The total dijet cross section, <r, can be factorized into 
the partonic cross section, a, and the parton distribution 
functions (PDF), P(x\^, Q 2 ), in the proton and antipro- 
ton: a = xiP(x\, Q 2 ) X2P(x2, Q 2 ) a. The partonic cm. 



energy, y/§, equals ■ s /x\X2S. For sufficiently large values 
of x\ and x%, any large a s ln(s/Q 2 ) terms in a correspond 
to large a s ln(s/Q 2 ), which are of the order of a s Arj, and 
factorize in a . Using the BFKL prescription to sum the 
leading logarithmic terms a s ln(s/Q 2 ) to all orders in a s , 
results in an exponential rise of a with A77 H : 
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where «bfkl is the BFKL intercept that governs the 
strength of the growth of the gluon distribution at small 
x. In the leading logarithmic approximation (LLA), 
«bfkl is given by [0: 



«BFKL — 1 



a s (Q) 12 In 2 



(4) 



The predicted rise of the partonic cross section with Ar] 
is difficult to observe experimentally due to the depen- 
dence of the total cross section on the PDF. To overcome 
this difficulty, we measure the cross section at two cm. 
energies, ^/sa = 1800 GeV and <Jsb — 630 GeV, and 
take their ratio for the same values of x\, X2 and Q 2 . 
This eliminates the dependence on the PDF, and reduces 
the ratio to that of the partonic cross sections. The latter 
is purely a function of the A77 values: 
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Thus, varying ^/s, while keeping x\ % X2 and Q 2 fixed, 
is equivalent to varying Ar], which directly probes the 
BFKL dynamics. In addition, measurement of the ratio 
leads to cancellation of certain experimental uncertain- 
ties, and enables an experimental extraction of q;bfkl • 

In the D0 || detector, jets are identified using the 
uranium/liquid-argon calorimeters. These cover the 
range of | r] |< 4.1, and are segmented into towers of 
A77 x A(j> = 0.1 x 0.1 (cf> is the azimuthal angle). 

The data samples for this analysis were collected dur- 
ing the 1995-1996 Tevatron Collider run. Events were 
selected online by a three-level trigger system culminat- 
ing in the software trigger requirement of a jet candidate 
with E T > 12 GeV. The trigger was 85% efficient for 
jets with Et = 20 GeV, and fully efficient for jets with 
Et > 30 GeV. The integrated luminosity of the trigger 
was 0.7 nb _1 for the ^fs = 1800 GeV sample, and 31.8 
nb" 1 for the </s = 630 GeV sample §. 

Jets were reconstructed offline using an iterative fixed- 
cone algorithm with a cone radius of 1Z = 0.7 in (77, (j)) 
space ||. The pseudorapidity of each jet was correct- 
ed for small reconstruction and jet algorithm biases. 
The transverse energy of each jet was corrected in three 
stages: (i) Energy originating from spectator parton in- 
teractions, additional interactions, noise from uranium 
decay, and residual energy from previous pp interactions 
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was subtracted on average from the measured jet ener- 
gy (ii) The jet energy was corrected for the hadronic 
response of the calorimeter M ; (iii) The fraction of the 
particle energy that showered outside of the jet reconstru- 
ction cone was recovered, and the fraction of the energy 
reconstructed within the cone that did not belong to the 
original particle was subtracted |l(| . The average correc- 
tion for jets of E T = 20 GeV and M = 2.5 is (22.8 ±4.8)% 
at \/s= 1800 GeV; for jets of the same Et and \q\= 1.2 
the correction is (14.5 ± 4.0)% at 630 GeV. 

The event vertex was required to lie within 50 cm of 
the detector center; 93%(86)% of the events at 1800(630) 
GeV satisfied this requirement. To remove cosmic ray 
background, the imbalance in the transverse momentum 
of the event was required to be less than 70% of the lead- 
ing jet Et] more than 98% of the events at each cm. 
energy satisfied this requirement. To ensure good jet 
reconstruction efficiency and jet energy calibration, jets 
were selected with Et > 20 GeV and |?7|< 3. Back- 
grounds from isolated noisy calorimeter cells, accelerator 
beam losses, and electromagnetic clusters that mimic jets 
were eliminated by applying a series of jet quality crite- 
ria; 97% of the jets survived this final selection. 

The selected jets of each event were ordered in pseudo- 
rapidity. A minimum pseudorapidity interval of Ar/ > 2 
was required between the most forward and most back- 
ward jet. In the final samples, the most forward and most 
backward jets were found to have approximately the same 
Et- The values of x\, X2 and Q 2 were calculated from 
Eqs. (0) and @. Most of the data at y/s = 1800 GeV 
are within 0.01 < 11,2 < 0.30, and at 630 GeV, within 
0.03 < x\ t 2 < 0.60. The region of maximum overlap, 
0.06 < 2:1,2 < 0.30, was divided into six equal bins of x% 
and X2- Due to limited statistics, only one bin in Q 2 was 
used: 400 < Q 2 < 1000 GeV 2 . The dijet cross section, 
corrected for trigger, event and jet selection inefficiencies, 
was computed in each (xi,X2,Q 2 ) bin. 

The dijet cross section at low (2:1,2:2) is affected by the 
acceptance of the Et > 20 GeV and A77 > 2 require- 
ments. To avoid this bias, we require x% ■ X2 > 0.01. 
Similarly, the cross section at high (£1,212) is biased by 
the \r]\< 3 requirement, so that we require 0:1,2 < 0.22. 
A total often (xi,X2) bins satisfy both requirements. 

Multiple pp interactions during the same beam cross- 
ing, which, in principle, could distort the topology of 
the event and bias the cross section, were infrequent for 
the low instantaneous luminosity (£ < 10 30 (2 x 10 30 ) 
cm _2 s _1 at y/s — 1800(630) GeV) data used in this anal- 
ysis. Nevertheless, any possible luminosity effects on the 
dijet cross section were evaluated by measuring the cross 
section at 630 GeV from lower- and higher-luminosity 
subsamples. No significant difference was observed be- 
tween the two measurements. 

The dijet cross section is distorted by jet energy reso- 
lution. The resolution was measured as a function of jet 
pseudorapidity and Et, by balancing Et in events with 



only two jets back-to-back in (f>. For jets of Et = 20 GeV, 
the fractional E T resolution is 27%(14%) at \tj\ = 1.2(2.5) 
and y/s = 1800 GeV (measured from the jet data col- 
lected during the 1994-1995 92 pb _1 Tevatron run). At 
= 630 GeV, limited statistics prohibited the measure- 
ment of the resolutions in the whole Et and rj spectrum. 
In the regions where the measurement was possible, the 
resolutions at 630 GeV were found to be smaller than the 
resolutions at 1800 GeV by -1%. 

The distortion of the cross section was corrected us- 



ing the herwig 1 11 1 Monte Carlo (MC) event gener- 
ator, convoluted with the CTEQ4M Q PDF. In the 
MC events, the jet transverse energies were smeared us- 
ing the resolutions extracted from the 1800 GeV data. 
The E T > 20 GeV, 1 77 1 < 3 and A77 > 2 requirements 
were applied separately to the original fully-fragmented 
(particle-level) jets and to the -&r-smeared jets. Particle- 
level and smeared dijet cross sections were calculated in 
the same (xi,X2,Q 2 ) bins as in the data. Apart from nor- 
malization differences, the smeared herwig cross section 
at both cm. energies exhibits the same dependence on 
2:1,2 as the data. The ratio of the particle-level to the 
smeared MC cross section in each bin was used as an un- 
smearing factor to correct the data cross section for the 
jet energy resolution effects. The unsmearing correction 
for the dijet cross section is typically of the order of 10% 
at both cm. energies, whereas the unsmearing correc- 
tion for the ratio of the cross sections amounts to only 
6%. The difference between the measured resolutions at 
the two cm. energies was accounted for in the systemat- 
ic uncertainties. The unsmearing method was verified by 
using a smeared MC sample generated with ISAJET [|l3| , 
and comparing the isajet particle-level cross section to 
that obtained using our unsmearing procedure based on 

HERWIG. 

The dijet cross sections for Arj > 2 at yfs — 1800 and 
630 GeV in the selected (xi,X2) bins are shown in Ta- 
ble |. In each bin, the average values of Xi, 2:2 and Q 2 
are in good agreement, within the precision of our mea- 
surement, between the two cm. energies. This ensures 
the cancellation of the PDF in the ratio of the cross sec- 
tions. Also shown in the Table are the values for the 
BFKL intercept, q^bfkl , extracted from the cross sec- 
tions and the average pseudorapidity intervals at 1800 
and 630 GeV in each (xi,X2) bin, using Eq. (||). 

The mean value of the ratios of the cross sections in the 
ten bins is equal to (R) = (tTi8oo/c63o) = 2.8 ± 0.3 (stat). 
The mean value of «bfkl is equal to 1.65 ± 0.05 (stat). 
The mean pseudorapidity interval, (A77), in the selected 
bins is equal to 4.6 units at 1800 GeV and 2.4 units at 
630 GeV. 

The largest sources of systematic uncertainties on the 
ratio of the cross sections and the BFKL intercept are 
the jet energy scale (yielding an 8% uncertainty on the 
ratio and 2% on the intercept) and the jet energy reso- 
lutions (7% on the ratio and 2% on the intercept). The 
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TABLE I. The dijet cross sections for A77 > 2 at y/s = 1800 
and 630 GeV and the extracted value of the BFKL intercept 
in each of the ten (xi,X2) bins. The minimum jet Et is 20 
GeV. The uncertainties are statistical. 
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individual components of these were evaluated for corre- 
lations between the two data samples. Additional sources 
of systematic uncertainties on the ratio and the intercept 
include the choice of the input PDF in the Monte Car- 
lo used for unsmearing (1% on the ratio, negligible on 
the intercept) and the uncertainty in the normalization 
of the luminosity (2% on the ratio and 1% on the inter- 
cept). The total systematic uncertainty amounts to 11% 
on the ratio of the cross sections and 3% on the BFKL 
intercept, yielding the final results: 
(R) = 2.8 ± 0.3 (stat) ± 0.3 (sys) = 2.8 ± 0.4, 
(o^bfkl) = 1.65 ± 0.05 (stat) ± 0.05 (sys) = 1.65 ± 0.07. 
Hence, for the same values of x%, X2 and Q 2 , the dijet 
cross section at large A 77 increases by almost a factor of 
three between the two cm. energies, corresponding to the 
increase of (A77) from 2.4 to 4.6 units. 

Several theoretical predictions can be compared to our 
measurement. Leading order QCD predicts the ratio of 
the cross sections to fall asymptotically toward unity with 
increasing Ar). For the Ar) values relevant to this analy- 
sis, the predicted ratio is i?LO = l-2 |l4| ]. 

The herwig MC provides an alternative prediction. 
It calculates the exact 2-^2 subprocess, including ini- 
tial and final state radiation and angular ordering of the 
emitted partons. Using the same (xi,X2,Q 2 ) bins as in 
the data yields -Rherwig = 1-6 ± 0.1 (stat). 

The LLA BFKL intercept according to Eq. (||) for 
a s (20GeV) = 0.17 (H) is a BFK L, lla = 1-45. For 
Ar/isoo = 4.6 and Ar te30 = 2-4, Eq. (§) yields 
^bfkl. lla = 1-9. It should be noted, however, that 



the leading log approximation may be too simplistic, and 
that exact quantitative predictions including the next-to- 
leading logarithmic jl5| corrections to the BFKL kernel 
are not as yet available. 

It is evident that the growth of the dijet cross section 
with Ar) (from (Ar)) — 2.4 to 4.6) is stronger in the data 
than in the theoretical models we considered. The mea- 



sured ratio is higher by 4 standard deviations than the 
LO prediction, 3 deviations than the herwig prediction, 
and 2.3 deviations than the LLA BFKL prediction. 

It should be noted that the x\ t 2 definitions of Eq. ([!]) 
have been kept the same in the data and in the theoretical 
calculations. Modifying these definitions to account for 
all jets in the event changes the ratio of the cross sections 
by less than 10%. 

Finally, the Ar) > 2 requirement was changed to 
Ar) > 1, and the analysis was repeated. For Ar) > 1, 
Eq. ([!]) yields x% ■ X2 > 0.005, which results in a selection 
of fifteen unbiased (x±,X2) bins. The mean pseudorapidi- 
ty interval in the selected bins is equal to 4.2 at 1800 GeV 
and 1.9 at 630 GeV. The average ratio of the 1800 and 630 
GeV cross sections in the selected bins was measured to 
be 1.8±0.1 (stat)±0.1 (uncorrelatedsys). The results are 
shown in Fig. |l| as a function of the mean pseudorapidity 
interval at y/s = 630 GeV. In the case of the Ar) > 1 
requirement, the observed ratio is once again larger than 
the exact LO and herwig predictions. It is interesting, 
however, that herwig exhibits the same qualitative be- 
havior as the data in that the ratio of cross sections de- 
creases as the Ar) requirement is relaxed, whereas the 
exact LO calculation predicts a very different trend. (A 
BFKL prediction is not shown for the case of Ar/ > 1 
since the pseudorapidity interval is not sufficiently large 
for the formalism to be meaningful.) 
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FIG. 1. The ratio of the dijet cross sections at y/s = 1800 
and 630 GeV for A77 > 1 and A77 > 2. The minimum jet Et 
is 20 GeV. The inner error bars on the data points represent 
statistical uncertainties; the outer bars represent statistical 
and uncorrelated systematic uncertainties added in quadra- 
ture. The error bars on the herwig predictions represent 
statistical uncertainties. The LO and BFKL predictions are 
analytical calculations. 

In conclusion, we have measured the dijet cross section 
for large pseudorapidity intervals at y/s = 1800 and 630 
GeV, and the ratio of the cross sections for the same 
values of x\, X2 and Q 2 at the two energies. The latter 
corresponds to the ratio of the partonic cross sections 



6 



for different values of Arj. The measured partonic cross 
section increases strongly with Ar), more strongly than 
expected on the basis of any current prediction. 
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